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We present a class of exact ground states of a three-dimensional periodic Anderson model at 
3/4 filling. Hopping and hybridization of d and / electrons extend over the unit cell of a general 
Bravais lattice. Employing novel composite operators combined with 55 matching conditions the 
Hamiltonian is cast into positive semidefinite form. A product wave function in position space allows 
one to identify stability regions of an insulating and a conducting ground state. The metallic phase 
is a non-Fermi liquid with one dispersing and one flat band. 



The periodic Anderson model (PAM) is the basic mi- 
croscopic model for the investigation of heavy fermion 
and intermediate valence systems, e.g., compounds con- 
taining elements with incompletely filled / shells such 
as Cerium or Uranium [1]. In its simplest form the 
PAM describes strongly correlated, dispersionless / elec- 
trons which couple, via a local hybridization, to non- 
interacting conduction (d) electrons hopping between 
nearest neighbor sites. For real systems this is certainly 
an oversimplification since there is experimental evidence 
for (i) a weak, but finite, dispersion of the / electrons (es- 
pecially in Uranium compounds) [2], (ii) non-local con- 
tributions to the hybridization, and (iii) hopping of the 
d electrons beyond nearest neighbors [3]. 

Recently, investigations into the origin of the dramatic 
volume collapse at the a — > 7 transition in Cerium 
have drawn attention to the possibility of a Mott metal- 
insulator transition in the PAM. [4-6] Although such a 
transition is usually associated with the half-filled Hub- 
bard model with nearest neighbor hopping, there exist 
remarkable similarities between the two models [7] , espe- 
cially if both the d electron hopping and d — f hybridiza- 
tion in the PAM connect nearest neighbor sites. [7-10] 
These results show that the spatial range of the hopping 
and hybridization in the PAM are very important, even 
on a qualitative level. 

In this situation exact results on the existence of in- 
sulating and metallic phases in the PAM and their de- 
pendence on a general set of hopping, hybridization and 
interaction parameters are particularly desirable. Exact 
results for the PAM are very rare since, in contrast to the 
Hubbard model, there does not even exist an exact solu- 
tion of the PAM in dimension D — 1. On the other hand 
it has been possible to construct exact ground states of 
the PAM in certain regions of parameter space, namely 
for infinite repulsion of the / electrons [11-13], and for 
finite repulsion in low dimensions {D — 1,2) [14,15]. 

In this Letter we show that it is possible to construct 
exact ground state wave functions describing metallic and 



insulating phases of the PAM at non-integer electron fill- 
ing even in dimension £) = 3. In particular, we explicitly 
demonstrate (i) the insulating and conducting nature of 
the solutions, (ii) the presence of strong variations in the 
compressibility of the system when leaving the insulat- 
ing phase, and (iii) the non-Fermi liquid nature of the 
metallic phase. 

We consider a general Bravais lattice in _D = 3 with 
a unit cell / defined by the primitive vectors {xt-}, r — 
1,2,3. The Hamiltonian of the PAM has the form H = 
Ho + UD^ where = describes the local 

Coulomb repulsion between the / electrons (C/ > 0). The 
one-particle part may, in general, be written as 

i,(7 ^ r 

+fLdi+r,a) + H.c] + Uo(4/i. + H.c.) + Efhl I , (1) 

where the terms with f^'^ represent the kinetic energy 
of d and / electrons due to hopping between two sites i 
and i -I- r, 14 is the hybridization of d and / electrons at 
sites i and i -I- r, Vb is the local hybridization, and Ef \s 
the local / electron energy; here r 7^ 0. The amplitudes 
f^'f are real, but V^, Vq can, in principle, be complex. 
In our investigation the hopping and hybridization of the 
d and / electrons extend over the unit cell of a general 
Bravais lattice (Fig. 1). To avoid multiple counting of 
contributions to (1) by the H.c. term, the vector r must 
be properly defined. To this end the sites within /i, the 
unit cell defined at site i, are denoted by r/j = i -I- r^^^, 
with Vap^ = axi -1-/3x2 -1-7x3; a, /?, 7 = 0, 1. As shown in 
Fig. 1 the eight sites r/j can be numbered by the indices 
ri(a,/3,7) — l-|-a-|-3/3 + 47 — 2a/3 without reference 
to Ii. In the following we use orthogonal x,- vectors for 
simplicity. Then r = Ya'p'-,' — Vajs^, with n{af, l3f,jf) > 
n{a, (3,"f), connects any two sites within a unit cell [16]. 

We now introduce a superposition of operators creating 
b {— d or /) electrons inside every unit cell Ii as 
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1+X1+X2+X3 
n=7 
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n=l -"Xi jj^2 

FIG. 1. Unit cell (/) of an orthorhombic lattice at an ar- 
bitrary site i showing the primitive vectors Xr and indices 

n of the sites in I. Arrows depict some of the hopping and 
hybridization matrix elements (J = t, V) extending over I. 
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b=d,f a,/3,7=0 

E Kbbl + <bbU.,,a + - + «8,6&U,+X3, J. (2) 

b=dj 



with ^ for all n. Although {A\^,A\,^,} = 

{Ai^,Ar^,} = and {Aj^, Alj^Kd + K f, where Kb = 
Xln=i kn.&P = "-(c^j/JjT), the composite opera- 

tor A\^ does not obey canonical anti-commutation rules 
since {Aja,A\,^,} ^ for / 7^ /'. Due to the transla- 
tional symmetry of the lattice the prefactors a* are the 

same in every unit cell. Making use of this fact, H can 
be cast into the form 



(3) 



where Pi = n(^n(^ - -n(^ + l, Eg= K^N + UNa - 



2N\(2Kd — Ef), and N and is the number of electrons 
and lattice sites, respectively. For (3) to reproduce (1) 
the prefactors o* in A\,^ must be expressed in terms of 
the microscopic parameters i^, t(, Vr, V*, Vb, Vq , Ef, 
U, for all r G /i, taking into account periodic boundary 
conditions. This leads to 55 coupled, non-linear matching 
conditions [17] which can be written in compact notation, 
with b, h' = d, f, as [18] 
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/3i,/32,/33=-l 



(11 ^ft , a J<+,ba« 



nb,b' 



(4) 



Apart from the constant term Eg in (3) H is positive 
semidefinitc operator. A state |\l/g) fulfilling the condi- 
tions Pi|*g) = and A\,J^g) = for all i will then 
be an exact ground state of H with energy Eg. We note 
that Pi assumes its lowest eigenvalue, 0, when there is 
at least one / electron on site i. Therefore a state of 
the form |*g) ~ Hilfi-^^lO) > where the operator = 
A*iT/iT (with arbitrary coefficients fii^) creates one 



/ electron on every site i, fulfills the first condition. Fur- 
ther, a state of the form [*<,) ~ n^'\(^liT^/ii)l^) ^"^^^^^ 
the second condition. Consequently, the (unnormalized) 
product state 



\^g) = n [AiAa^i] io) 



(5) 



i=l 



has at least one / electron on every site i, the prod- 
uct riila^/ T^/ i creating at most two more (d or / ) 
electrons on i. Clearly, has the desired property 

H\^g) = Egl^g) 8,11(1 IS thus au exact ground state of H 
with energy Eg . Since the product of the three operators 
in I'^g) creates N = 3Na electrons, the ground state is 
3/4 filled, i.e., there are on average three electrons per 
site. The arbitrariness of /ii ^ implies a large spin degen- 
eracy of |\E'g) which is globally paramagnetic. Neglecting 
the trivial 25-1-1 multiplicity of the spin orientation, the 
ground state is 7VA/2-fold degenerate [19]; all degenera- 
cies are still contained in (5). We note that is only 
a ground state of the interacting system since the opera- 
tor Fi in l^fg) is required to enforce t^X^i Aj^'g) = for 
U > 0. Hence |\E'g) is not connected in any simple (e.g., 
perturbative) way to the non-interacting {U = 0) ground 
state of the model. 

The physical nature of |\I'g) depends on the values of 
the coefficients a* ^ in (2) which are solutions of (4) for 
given microscopic parameters. We will now identify local- 
ized and itinerant ground states, and discuss their phys- 
ical properties. 

1. Localized ground state. The coefficients a„ j, may be 
chosen in such a way that \'9g) has exactly three electrons 
per site. To see how this can be achieved we take a look 
at a typical factor j^Fjt, entering in 

where j" is a common site of the two unit cells /j, ly. 
Since Fy, always creates one / electron on j", the product 
of the four unit cell operators should create only two more 
electrons (one f and one |) on j". Therefore terms of the 

form 4'T/j"T4'i/j"i(<i,d<2,/ -<i,/<2,d)' which are 
also generated and which lead to more than two addi- 
tional electrons on this site must be prohibited. This can 
be achieved by choosing a* ^/ a* j = p for all n. It follows 
from (4) that along the diagonal of a unit cell with end 
points n,n' the relation a'^jan',d = c-n^^n'j holds, im- 
plying p and therefore also the hybridization amplitudes 
to be real. A solution of (4) of this form is obtained 
for t^ = t\,, i.e., for equal hopping along x,y,z. For the 
coefficients chosen in this way the ground state becomes 

I*;-) = HE l^i^Pipdlidhfl + //t/u4)] |0>- (6) 

i CT 

Denoting ground state expectation values by (...) the lo- 
cal /, (i-elcctron occupations are found as (n/) = (1 -|- 
2^)/(l + z), {hf) = (2 + z)/(l + z), where z = \t(/ti\ is 
a measure of the nearest-neighbor hopping amplitude of 
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the / electrons. Since the latter can be expected to be 
much smaller than that of the d electrons {z « 1), the 
/ (d) electron occupation per site is found to be close to 
one (two). Hence there exist local moments on most of 
the / sites. We see that (fii) = {n() + (nf) = 3 for all 
i, i.e., the electron distribution is indeed uniform. The 
localization is due to a subtle quantum mechanical inter- 
ference between the hopping and hybridization processes 
of the electrons. Therefore the nature of this localized 
state is quite non-trivial. 
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FIG. 2. Ground-state energy of the localized state ex- 
pressed in terms of {Eg/NA + ?7)/ |tf | as a function of 
next-nearest neighbor hopping of d electrons, \ti/tf\, for dif- 
ferent values of / hopping z = \t{/tf\; here |Vi/Vb| > 1/2. 

By separating the Hamiltonian H into an itinerant 
part Hitin = Ho- Ei,a[^oidlfi<^ + ^-c.) + Efhl] = 
'^r ^itini^) ^nd a complementary localized part Hioc = 
H - Hiun- and using {b\J}'^„,) = for i ^ j and 

b,b' = d,f, one finds {Hitin(r)) = and (Hioc) = Eg. 
This clearly illustrates the localized, and thus insulating, 
nature of the ground state whose energy is obtained as 
Eg = Na[{1 - '2z)Kf/z - U] (Fig. 2). The parameter 
space in which the localized ground state |\E'joc) is stable 
is represented in terms of the variables Ef, U, t(/tf, tf 
by the surfaces Ji, I2 in Fig. 3. The stability region is 
seen to extend through the phase diagram from weak to 
strong interactions U . 

For y = \t2/tf \ > 1/2, i.e., rather large next-nearest 
neighbor hopping of d electrons, the localized state |\l/(oc) 
ceases to be the ground state. Apparently, at yc = Ifi 
a different, most probably itinerant, phase becomes sta- 
ble. We note that the ground state energy Eg{y) has a 
finite value at yc, but infinite slope (see Fig.2). This ex- 
act result has a direct physical interpretation. Namely, 
since the size of the hopping element may be tuned by 
pressure, the infinite slope of Eg at y = yc is expected 
to correspond to an anomaly in the compressibility at a 
critical pressure Pc- Such a feature is indeed observed in 
some heavy- fermion systems [21]. 

2. Itinerant ground state. The localized ground state 
discussed above has exactly three electrons per site. In 
general, the intersite hopping and hybridization will lead 



to a variable number of electrons per site. For such states 
(Hitin) 7^ 0, implying a non-zero conductivity of the sys- 
tem as expressed by the sum-rule divRear^ri^) = 

-Tre'^\xr\'^/{2nV)J2^.^,^^_^g{HiUnir)), where ar,T{t^) is 
defined by the Kubo formula [20]. An itinerant ground 
state is obtained, for example, by choosing p* = 
—Pm\Pn\ = IpI) corresponding to imaginary p and, 
hence, imaginary hybridization amplitudes T4. Whether 
Vr is real, complex or imaginary depends on the lin- 
ear combination of the corresponding electronic orbitals 
[7,15,22-24], and hence on the lattice symmetry. For ex- 
ample, Vr may be tuned from real to imaginary by in- 
troducing axial distortions of D^h symmetry to an un- 
derlying Oh lattice symmetry [23]. Therefore such a 
solution requires anisotropic hopping and hybridization 
amplitudes. The itinerant ground state discussed here 
emerges if the hybridization on the same site and in the 
basal (xy) plane vanish. The anisotropy in the hopping 
starts at the level of next-nearest neighbor amplitudes. 

To describe the itinerant case a k-type representation 
is more suitable. Denoting the Fourier transforms of Ai^„ 
and bia by and Sku) respectively, (2) takes the form 
/ '^kb^ka! expressions for the coefficients 
oj^f, will not be reproduced here. We can now define 
new canonical Fermi operators Cs^ucr, S = 1,2, where 
Ci^ka = R-l/^Aka, with = J2b l«kbP, and C2,ko- is 
determined by the anti-commutation rules between Cj^.ko- 
and k<7- It follows that - J2i,a A.aM'y + KdN = 

Ek,a i^d - ^k ^)C'I,k<7C'l,k<T + KdCl i^^C2,ka = Hg , 

such that (3) can be written as H = Hg + UP + Na[U - 
2{Kd - Ef)]. In the ground state, using P\^g) = 0, the 
Hamiltonian H therefore reduces to Hg. Thus we suc- 
ceeded in diagonalizing H for the ground state. There are 
two bands, the lower one having a dispersion — R^^, 
while the upper one is dispersionless ("flat"); the Fermi 
energy is Ep = K4. Such a band-structure around Ep 
has been observed in experiment [25]. 

The momentum distribution of the d, f electrons be- 
comes 



(E^LM = (2|ak6p + lak.b'TtbH-Rk, 



(7) 



with nk = J2b=d / '^k = 3. Since the coefficients ak6 are 
regular functions of k this also holds for n^, n{^, and rik. 
Consequently, the momentum distributions of the elec- 
trons in the interacting ground state has no discontinu- 
ities. Since the ground state is paramagnetic and metal- 
lic, the system is a non-Fermi liquid. This is a conse- 
quence of the degeneracy of electrons in the upper band. 
In terms of the C^^ka fermions one finds {Cl ^^^Ci^i^^r) = 1 

and (C"! kCT^2,kcr) = 1/2 (upper band half filled). 

The itinerant solution can be generalized to fill- 
ings beyond 3/4 by inserting the operator = 
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nj^iEr=i aj*(E6=d,/> ^bab]^)] into (5) next to |0); here 
ftji) Cbcr are numerical coefRcients. This operator intro- 
duces M < A/a additional particles into the ground state. 
It allows one to calculated the energy Eg for different par- 
ticle numbers. In particular, one finds /x+ = Eg{N + 2) — 
Eg{N + 1) = Kd, II- = EgiN + 1) - EgiN) = Kd, i.e., 
fi~^ — IJ,~ = 0. Therefore the system is a conductor. [26] 
The stability region of the conducting state corresponds 
to the surface C in Fig. 3. 




FIG. 3. Surfaces in parameter space representing stability 
regions of the insulating ground-state (/i for |Vi/Vb| < 1/2, 
I2 for |Vi/Vb| > 1/2), and the conducting ground-state (C). 

In summary, we derived the first exact ground state 

solution of a f/iree-dimensional periodic Anderson model 
with finite hopping and hybridization of d and / elec- 
trons in the unit cell. This was achieved by (i) cast- 
ing the Hamiltonian into a positive scmidefinite form us- 
ing composite operators in combination with 55 coupled, 
non- linear matching conditions, and (ii) constructing a 
product wave function of these composite operators in 
position space. For real hybridization amplitudes we ob- 
tained an insulating ground state whose compressibility 
diverges at the boundary of the stability region. By con- 
trast, for imaginary hybridization amplitudes we identi- 
fied a conducting non-Fermi liquid state consisting of one 
dispersing band and one (upper) flat band. The stabil- 
ity regions of the two ground states extend through an 
unexpectedly large region of parameter space, e.g., from 
weak to strong interactions U. 

By modifying the structure of the composite operators 
it is possible to vary the stability regions of these ground 
states, and also to describe magnetically ordered phases. 
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